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1 Introduction 



One ingredient of the contemporary Standard Model is QCD which is beheved to 
describe the strong interaction between quarks and gluons. The theory has six 
free mass parameters for the six quark species and one couphng parameter. After 
fixing these free parameters of the QCD Lagrangian at some reference energy 
scale, the theory in principle predicts presently known effects in which only the 
strong interaction is involved. Since QCD possesses the property of asymptotic 
freedom, perturbation theory is applicable in the high energy regime {E ^ IGeV), 
but for low energies lattice QCD, the only known systematic non-perturbative 
approach, is required. Here a certain number of physical quantities, for example 
meson masses, have to be identified with experimental inputs to determine the 
free parameters and in particular to gauge dimensionful quantities on the lattice 
in powers of GeV. Only then predictions become possible. As an example, the 
light hadron mass spectrum was studied by different collaborations in quenched 
as well as in full QCD with up to three dynamical fermion species [1, 2, 3, 4, 5, 6]. 
The theoretical predictions for the spectrum in full QCD are compatible with the 
experiments at the error-level achieved so far. 

At a first glance the perturbative and non-perturbative formulations seem to 
be disjoint and have their own domains of applicability. In a long-term project, 
the ALPHA collaboration has been developing methods and tools to connect these 
regimes by computing the parameters of the perturbative domain starting from 
the non-perturbative formulation [7]. In this context, the step scaling function 
which we will discuss later plays a key role. For different numbers of flavors, the 
step scaling function [8, 9, 10, 11, 12, 13] which can be understood as an integrated 
form of the QCD /3-function and the running coupling have been examined and 
determined successfully (see [14, 15, 16] for recent publications). In this paper, 
we want to calculate the step scaling function and the running coupling with four 
flavors of 0{a) improved Wilson quarks. For the present work, the improvement 
coefficient Csw which is essential for the 0{a) improvement with Wilson quarks was 
determined in a preceding paper [17]. We exploit this result in the simulations 
that we report in the following. 

The paper is organized as follows. In sections 2-5 we summarize some back- 
ground about the nonperturbative definition of a coupling constant in the Schro- 
dinger functional finite volume renormalization scheme and about the step scaling 
technique. In section 6 we present our new raw data and their analysis. We will 
arrive, for four flavors, at a value for the dimensionless combination ALmax- While 
A parameterizes the coupling at high energy, Lmax is an unambiguously defined 
length in the hadronic regime. Its value in GeV remains to be determined by large 
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volume simulations. 



2 The finite size strategy 

The fundamental parameters of QCD, i.e. the coupling and masses of the quarks, 
depend on a scale fi. To make contact with perturbative QCD we need to know 
their values for large /i in the domain of asymptotic freedom. If such a computa- 
tion is attempted on a single lattice where at the same time hadronic scales are 
measurable with small finite volume and discretization effects one has to satisfy 
to a reasonable precision the string of inequalities 

1111 , , 

^> ~ 77TTTTT7 > - ~ T777TT7 > «• (2-1) 

0.14 GeV lOGeV ^ ^ 

In total this implies L/a ^ 70 which is impossible to satisfy for some time. A 
way to circumvent these difficulties was hence proposed by Liischer et al. in 
[18]. The idea is to exploit the universality of the finite volume continuum limit 
to perform intermediate renormalizations with the finite size as a renormalization 
scale, /i = 1/L. In these steps only L ^ a has to be assured to render cutoff effects 
small. By such simulations one can nonperturbatively determine the change of 
the coupling constant in the continuum limit under scale changes from L to L/2. 
By repeating such steps we start at Lmax ~ C(0.5fm) and arrive after k steps at 
H = With a sufficient number of steps this will be in the perturbative 

regime where we can make contact with the A parameter of an arbitrary scheme 
and with the scheme independent renormalization group invariant quark masses by 
applying perturbative formulas. If one finally succeeds in computing the precise 
value of //max in units of a mass or decay constant, then all references to the 
intermediate finite volumes will have disappeared from the final result. More 
details may be found in [7]. 



3 Running coupling 

As mentioned before, the coupling and the quark masses are renormalization scale 
dependent and run with energy. Therefore the quoted values of the quark masses 
and the world average of as, both in the MS scheme, in the Particle Physics 
Booklet [19] are referred to a particular reference scale (// ~ 2 GeV for the masses 
and /i = Mz for the coupling). 

From the theoretical point of view, the running of the QCD parameters is 
described by the renormalization group equation. However, a physical observable 



3 



G should have no reference to a particular renormalization scale yU. This fact is 
expressed by the Callan-Symanzik equation 

rrii 



dg 



d_ 

' dfrij 



0. 



(3.1) 



In words: for any change of fi there are accompanying modifications of g and rrii 
such that the physics is unchanged. The implied scale dependence of the coupling 
is given by /3{g) in (3.1) 

.te)^.^. (3.2) 

For high energies or for weak couplings, the /3-function has the following asymp- 
totic expansion 

9^0 



(3.3) 



The first two (1- and 2-loop) coefficients bo and bi in (3.3) are scheme independent 
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(3.5) 



The 3-loop coefficient does depend on the scheme and in the Schrodinger functional 
scheme, which we will need later, 62 has been given in [24] 

62 = — ^ [0.483(7) - 0.275(5)Arf + 0.0361(5)iVf2 - 0mi75(l)Nf] . (3.6) 

[4:71 j'^ 

We restrict ourselves to mass-independent schemes where all renormalization con- 
ditions are imposed at vanishing quark masses. Examples are the MS scheme of 
dimensional regularization and the Schrodinger functional scheme used here. For 
Nf < 16, the /3-function (3.3) is negative at weak coupling and the integration 
of (3.2) results in a coupling which decreases with increasing energy. In other 
words, the quarks behave like free particles in the high energy regime (asymptotic 
freedom). But for iVf > 16, there is a sign change and the property of asymptotic 
freedom is lost. The relation of g{fi) to the A parameter of QCD is given by the 
following solution of the Callan-Symanzik equation 

[bofifJ')] "^exp ' 
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dx 
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(3.7) 



As is well-known, A is scheme- dependent but the transformation to other schemes 
follows from the relation between their couplings at one loop accuracy. 
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4 Coupling in the SF scheme 



The finite-volume scheme which we use for our simulations is the Schrodinger 
functional (SF) scheme. In the following, we will only give a brief reminder of the 
main features and properties of this scheme which are discussed in detail in many 
papers, for example [20, 21, 22, 23, 24]. 

The original proposal was made in [20]. The construction was guided by 
several design criteria. The coupling constant was to be defined independently of 
perturbation theory, but its evaluation had to be manageable both in lattice per- 
turbation theory (also beyond one loop) and by numerical simulation. Moreover, 
small lattice artefacts were demanded. These requirements are not easy to fulfill 
and have led to the Schrodinger functional which is the Euclidean propagation 
kernel of a field configuration C at time xq = to another field configuration C 
at time xq = T. The free energy or effective action F of such a system is given by 

Z[C",C] = exp{-F} = [ ^[[/,^,^]exp{-^[?7,^,^]}. (4.1) 

The action S is defined as in [17]. The choice of Ck and C[ of the boundary gauge 
fields 

U {x, k) 1^0=0 = exp{aCfc}, (4.2) 
f/(x,fc)U„=T = exp{aCa (4.3) 

is largely arbitrarily. After some optimization with regard to lattice artefacts 
constant Abelian fields turned out to be appropriate. We will also adopt this 
choice of boundary fields parameterized by the scale L and two dimensionless 
real parameters 77 and v [8]. Numerical simulations of the Schrodinger functional 
showed that the choice z/ = leads to small statistical errors for the coupling in 
simulations. Therefore, we will also set v to zero. 

The Schrodinger functional coupling is now defined by the response of F to 
a variation of the boundary fields around point 'A' of [8] via the parameter 77, 



dr] 



r)=0 
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(4.4) 



where A; is a normalization constant chosen such that the perturbative expansion 
of F' begins with the bare coupling at tree level [8]. It should be noticed that 
the only external scale which appears in the definition of the coupling is the box 
size L, i.e. recursive finite size techniques can be used for the investigation of the 
evolution of the coupling. 
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The key quantity V in the definition of the couphng is an observable which 
can be calculated easily through Monte Carlo simulations. Taking the derivative 
of (4.1) results in 

r'=-^ln|y" ^[?7,^,^]exp{-5[t/,^,^]}| (4.5) 

= 1| m^,^][^ + ^)eM-S[U,^,m (4.6) 



dScX /dS- 



dr] / \ drj 



(4.7) 



Explicit expressions for both expectation values in (4.7) can be found in [11, 25]. 
The calculation of the renormalized coupling involves expectation values of a 
local operator and no correlation functions. Therefore the numerical evaluation 
on a computer is straight forward once configurations are available. The relation 
to the QCD coupling in the Schrodinger functional scheme asF is given by 

asF(;u) = ^^, /i = l/L. (4.8) 

To complete the definition of the coupling, also the boundary conditions for 
the quark fields have to be specified. We follow exactly [23]. In particular the 
angle entering the spatial periodicity of the quarks is chosen as ^ = 7r/5 since this 
value is advantageous for the numerical simulations [11, 23]. 



5 The step scaling function 

The concept of the step scaling function which was introduced in [18] has proven 
to be a very useful recursive technique to scale the coupling to high energies. As 
we discussed before, in our finite-volume scheme, the energy scale /i is identified 
with L~^. Hence the renormalization group function /3 (3.2) describes how the 
coupling changes if the box size is varied infinitesimally. The step scaling function 
a{s,u), in comparison, gives then a description how the coupling behaves when 
the box size L is scaled by a factor s 

fisL)=ais,f{L)). (5.1) 

The function a{s, u) can be regarded as an integrated form of the renormalization 
group /3-function. With the help of the step scaling function (5.1) the coupling 
can be traced to scales 2~^Linax (small box sizes, high energies) starting with an 
initial value L = Lmax- The value s = 2 is commonly used [11, 18, 26] in the 
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application of a{s, u) and we will also make this choice (from now on s = 2 and 
we will drop the argument s). The relation between the renormalization group 
/3- function and the step scahng function a{g'^{L)) is given by 



-21n(2) = I 

J u 



For small values of the coupling u, the step scaling function has the following 
perturbative expansion 

a{u) = u + SqU^ + Siu^ + S2U^ + . . . (5.3) 

where the coefficients are given by 

So = 26oln(2), (5.4) 
si = [26oln(2)]' + 26iln(2), (5.5) 
S2 = [2bo \n{2)f + lOboh [\n{2)f + 262 ln(2). (5.6) 

Both eq. (5.2) with the truncated /3 function as well as the expansion (5.3) can be 
used as the perturbative approximation of the step scaling function. They differ 
from each other by higher order terms. We take the first option when we compare 
our Monte Carlo results with perturbation theory. 

To study (j{u) numerically one starts with choosing several lattice sizes L/a 
and tuning the bare coupling and the hopping parameter k in such a way that 
the value of the renormalized coupling g'^{L) reaches some chosen value g'^{L) = u 
and the quark mass vanishes since we are in a massless scheme. The next step 
is to take L/a — 2L/a and simulate at the same bare parameters. The obtained 
coupling g'^{2L) from the latter simulations is a lattice approximation Ti{u,a/L) 
of the continuum step scaling function a{u). An extrapolation to the continuum 
of the data points at the same coupling u but growing L/a then leads to one 
value of the continuum function a{u). The procedure is repeated until a suitable 
range of u is covered. An appropriate functional description of the continuum 
step scaling function can be given in the end in the form of a suitable fit-function 
that interpolates the data. 

The lattice approximation T,{u,a/L) of the step scaling function contains 
remnant lattice effects of order a. The reason is that beside our nonperturbative 
value for Cg^ there are boundary improvement coefficients for which only pertur- 
bative estimates are available. We use those to the known order [8, 23, 24] which 
ensures that 0{a) cutoff effects in the step scaling function appear only starting 
at three-loop order. 
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The details of cutoff effects also depend on how the condition of a massless 
scheme is exactly implemented at a finite lattice spacing. As in the Nf = 2 
computation, we define the massless point on the smaller of the pair of lattices 
entering the step scaling function. More precisely, the (unrenormalized) PCAC 
quark mass, 

Udo + do)fA{xo) + CAadodofp{xo) 

'"(^»' = 2M^) 

at Xq = "2", 

mi = m , (5.8) 

is set to zero. The definition of the correlation functions /a, /p is found e.g. in 
[17, 27]. The second reference also contains a discussion of the influence of choosing 
the massless point in a way which differs from (5.8) at finite lattice spacing. For 
the improvement coefficient ca in (5.7) we used the 1-loop result [28]. 

In addition to the various improvement terms in the action and ca in (5.7) 
we exploit our knowledge of 5i and 82 from the perturbative calculation of 

6(u, a L) = ' \ — = 6i(a L)u + 62(a L)u^ + . . . , (5.9) 

a{u) 



Si{a/L) = 6io{a/L) + 6u{a/L)N{ (5.10) 
62ia/L) = 52Q{a/L) + 52i{a/L)N, + 522{a/L)Nf (5.11) 

with the coefficients taken from [26]. For A'^f = 4 we have 



L/a 


5i 


52 


4 


-0.0102 


0.0073 


6 


-0.0045 


0.0013 


8 


-0.0024 


0.00013 



Using these data we form the lattice step scaling function [29] 

^ ' ' ' l + 5i{a/L)u + 52{a/L)u^ ^ ' 

which we expect to have smaller overall cutoff effects. They still start at order 
a X but terms of order a™ x vP" are removed for all m and for n < 3 (in fact 
non-perturbatively in a). As mentioned previously, the order a x -u^ terms are due 
to the only perturbatively known boundary improvement terms. Their influence 
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was explicitly checked for iVf = 2 and found to be minor [11], such that also here 
we assume that the step scaling function converges effectively at a rate 

S(2)(n,a/L) =a(n) + C(a2). (5.13) 

Once the non-perturbatively determined continuum step scaling function <t{u) 
is known, the running of the coupling can be computed easily. For this purpose, 
one solves the recurrence 

Ui = Cr(Ui+i), Z = 0, . . . , n, Uo = timax = (-^^max) , (5.14) 

where Wmax is chosen such that the corresponding scale L^ax is in the hadronic 
regime of QCD. Proceeding in this way the coupling can be obtained over a wide 
range of energies. At a sufficiently large energy /i = 2'^/Lmax {k < n), perturbation 
theory can be applied for determining the quantity ALmax using (3.7) with the (3 
function truncated at 2-loop and 3-loop respectively. 

6 Numerical computation and results 

6.1 Simulation parameters and raw data 

The choice of the improvement coefficients q , Ct and Cgw were as discussed in our 
preceding paper [17]. The range of /3 for our simulations was limited by the validity 
range of Csw with four flavors (/3 < 5.0) [17]. Since our code allows only an even 
number of lattice points in each direction and lattices beyond L/a = 16 are too 
time-consuming for our present resources, we chose lattices {T = L) L/a = 4,6,8. 
We picked a sufficient number of values of /3 e [5.0, 9.5] to adequately map out a 
range g'^{L) = 0.9 . . . 2.7, g^{2L) = 1.0 .. . 3.5. With some tuning of the hopping 
parameter n the quark mass was kept small enough (|miL| < 0.005) such that 
mass-effects in the step scaling function are negligible. 

We performed our simulations on 4-5 crates of the apeNEXT machine at 
DESY Zeuthen over a period of about a year. The raw data are listed in the 
appendix. 

6.2 Analysis of data 

The computation of the step scaling function on the lattice as described in section 
5 requires g^{L) to be fixed to certain values u while the resolution a/ Lis changed. 
Previously this was realized by tuning (3 for each pair u,L/a [8, 9, 11]. Instead 
we here follow the more convenient proposal of [13] to pick a sufficient range and 
number of bare couplings for each considered L/a and interpolate the running 
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Figure 1: Data points for L/a = 8 and L/a = 16 and their interpolations. The 
parameter n in (6.1) is set to three. The hopping parameter k, was tuned only on 
the small lattices L/a = 4,6,8 such that the PCAC mass (5.8) vanished. 

coupling g'^{P,L/a) with a smooth function of /3. Afterwards the function allows 
access to any value of in the covered range. As an interpolation we took 

-1 

(6.1) 

motivated by perturbation theory. Note, however, that we do not fix the known 
perturbative expansion coefficients; we do not even require co,L/a = 1- Somewhat 
different forms and values n were checked and we verified that our results do not 
depend on the details of these interpolations. An example is shown in figure 1 for 
a pair of lattice sizes, namely L/a = 8 and L/a = 16. The coefficients Cm,L/a were 
determined by a standard fit. From the interpolated couplings we computed 
T,{u,a/L) for a number of values u starting from an initial coupling ^initial = 0.9 
and following roughly a sequence given by the recursion (5.14). 
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1.1787 


4 


1.286(2) 


1.287(2) 


2.2003 


4 


2.650(11) 


2.617(11) 




6 


1.297(5) 


1.301(5) 




6 


2.688(17) 


2.698(17) 




8 


1.298(6) 


1.301(6) 




8 


2.684(24) 


2.697(24) 


1.2972 


4 


1.428(3) 


1.430(3) 


2.6870 


4 


3.462(22) 


3.378(21) 




6 


1.442(6) 


1.448(6) 




6 


3.507(40) 


3.517(40) 




8 


1.446(7) 


1.450(7) 




8 


3.477(44) 


3.496(44) 



Table 1: Results for S(m, a/L) and S^^''(m, a/L) for different lattices and couplings 
u. 



6. 2. 1 Error propagation 

The uncertainties of the initial MC data which can be found in the appendix 
are statistically uncorrelated. For the purpose of propagating their errors, let us 
collect them in a vector x with components Xi. The above fit function g'^{P,L/a) 
may then be regarded a function f{x) of the initial data^. The error 6f of / is 
simply estimated by 

Since / is a relatively complicated function, it is convenient to avoid computing 
the derivatives df /dxi analytically. Instead we just compute them by a symmetric 
finite difference, with the step in Xj given by 5xi itself. This convenient method 
is applied for estimating the errors of all quantities derived from our data in 
the following. If desired also the correlation matrix of the errors of different 
observables can be obtained this way. 

'^We neglect that in the way we determine the interpolation of there is also a dependence 
on the uncertainties (5a;,;. In the fit we could also replace the errors 5xi by a smooth predefined 
function of /3. The results would not change significantly. 
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6.3 Results 



Our result for both step scaling functions Il{u,a/L) and the perturbatively cor- 
rected S(2)(M,a/L) are listed in table 1. As one can see in the visualization of our 
data, figure 2, the cut-off effects seem to be very small except for L/a = 4. As a 
precaution against higher order cutoff effects, we thus excluded the data set of our 
coarsest lattice from our analysis leading to the continuum step scaling function. 
We carried out three different analysis. 

• Constant fit: A fit of T,^'^\u,a/L) for L/a = 6,8 to a constant, for each u. 

• Global fit: A fit 

(u, a/L) = a{u) + pu^ {a/Lf. (6.3) 

with a separate, independent parameter cr(n) for each value u but a common 
parameter p modelling the cutoff-effects. 

• L/a = 8 data: Using directly a{u) = 1/8). 

The three different ansatze yield results which are in complete agreement with 
each other as seen in Table 2. The value of p in (6.3) comes out to be p = 
0.007(85) which is a good indication that cutoff effects are negligible in the data 
for L/a = 6, 8. 

In figure 2 we depict the continuum extrapolation with the constant fit. How- 
ever, as our final results we take just the L/a = 8 data. This is more conserva- 
tive and we are confident that the statistical errors dominate over residual cutoff 
effects. In particular these data agree with the L/a = 6 data and also using 
E(m, a/L) instead of T,^'^\u, a/L) has a negligible effect. 

Using a polynomial of degree five in u, we performed a constrained interpo- 
lation of the data in the fourth column in table 2. The coefficients up to were 
fixed by perturbation theory. Our fit 

a{u) = u + SqU^ + siu^ + 0.0036 - 0.0005 ^^^ < m < 2.7. (6.4) 

is shown in figure 3 (thick line). The perturbative step scaling functions are 
close to the one sigma range of the non-perturbative data points over the whole 
interval of the coupling u. We note also a peculiarity. Perturbation theory at 
3-loop lies below the 2-loop truncation of the /3-function and further away from 
the non-perturbative result. However, the 3-loop coefficient 62 (eq. (3.6)) in the 
Schrodinger functional scheme changes its sign between A^f = 2 and Af = 3 and 
it is rather small for A^f = 4. It is hence not unlikely that the 4-loop term would 
move the perturbative curve closer again. 
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Figure 2: Continuum extrapolation of the step scaling function (constant fit). 
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u 


constant fit 


(r[u) 
global fit 


L/a = 8 data 


0.9300 


1.002 (3) 


1.002 (3) 


0.997 (5) 


1.0000 


1.084 (3) 


1.084 (3) 


1.081 (4) 


1.0813 


1.182 (3) 


1.182 (4) 


1.181 (5) 


1.1787 


1.301 (4) 


1.301 (5) 


1.301 (6) 


1.2972 


1.448 (5) 


1.448 (7) 


1.450 (7) 


1.4435 


1.634 (5) 


1.634(10) 


1.637 (8) 


1.6285 


1.877 (7) 


1.877(16) 


1.880(11) 


1.8700 


2.209(10) 


2.207(27) 


2.212(17) 


2.2003 


2.698(14) 


2.694(49) 


2.697(24) 


2.6870 


3.507(30) 


3.50 (10) 


3.496(44) 



Table 2: Results of different fit procedures (as described in the text) for the 
continuum extrapolation of the step scaling function. 

Using the parameterization (6.4) of the step scaling function, we calculated 
the combination In(ALinax) starting from the highest coupling u^aax = ^?^(-^max) 
which was covered by our non-perturbative step scaling function and solved the 
recursion step (5.14) numerically to obtain the couplings Ui which correspond to 
the energy scales /i = 2YLmax- With the help of (3.7) and using the truncated 3- 
loop (3 function, we computed the values for In(ALinax) recorded in table 3. From 
the L/a = 8 results we quote 

In(AL^ax) = -2.294(83) at n^^x = 3.45 (6.5) 

as our final result. This determination of the A-parameter in units of Lmax has 
a precision of ~ 8%. It remains to gauge Lmax in physical units through a large 
volume computation. Therefore we here show the running of the coupling in the 
Schrodinger functional scheme in units of A. Figure 4 displays asF computed from 
the sixth column of table 3. We observe that upon the iterative application of the 
step scaling function the difference between the perturbative (using (6.5)) and the 
non-perturbative coupling is around a 3-sigma effect at the strongest coupling. 

7 Conclusions 

We computed the step scaling function of the QCD coupling in the Schrodinger 
functional scheme with four massless flavors. We used 0{a) improved Wilson 
quarks after first determining the Sheikholeslami Wohlert coefficient Csw [17]. The 
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Figure 3: The step scaling function for A'f = 2,4 and the perturbative results. 
The thickest line is the fit of our data points. The upper solid line is the N{ = 2 
result and the lower lines show the perturbative results. 

resulting cutoff effects are very small (figure 2, table 1) allowing for a continuum 
extrapolation. While the data are compatible with a constant for L/a > 6, we 
assume this form only for L/a > 8; the smaller lattices thus only enter the analysis 
by demonstrating that cutoff effects are small. We emphasize that this statement 
refers to the present level of statistical errors. If in the future statistical errors 
are further reduced, larger L/a will be necessary at the same time. It will be very 
interesting to see also the efficiency of computations with different regularizations 
of the Schrodinger functional as well as the corresponding test of the universality of 
the continuum limit. Most notably there are chirally rotated boundary conditions 
for the quarks [30, 31, 32] and staggered quarks [33, 34] for which results are 
expected soon. 

Already now, we do observe a small but significant deviation from 3-loop 
perturbation theory at the largest coupling reached in figure 4. It is about 10% 
(three standard deviations) and the Schrodinger functional coupling has a value 
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constant fit 


global fit 


L/a = 


8 data 


i 


Ui 


In(ALmax) 


Ui 


In(ALmax) 


Ui 


In(ALmax) 





3.45 


-2.028 


3.45 


-2.028 


3.45 


-2.028 


1 


2.660(14) 


-2.074(17) 


2.666(46) 


-2.066(56) 


2.660 (21) 


-2.073(26) 


2 


2.173(13) 


-2.117(24) 


2.179(45) 


-2.105(83) 


2.173 (20) 


-2.116(37) 


3 


1.842(11) 


-2.155(28) 


1.847(37) 


-2.141(97) 


1.842 (17) 


-2.153(44) 


4 


1.6013(90) 


-2.188(32) 


1.606(30) 


-2.17 (10) 


1.602 (14) 


-2.185(50) 


5 


1.4187(78) 


-2.217(35) 


1.422(25) 


-2.20 (11) 


1.419 (13) 


-2.213(56) 


6 


1.2748(70) 


-2.241(39) 


1.278(20) 


-2.23 (11) 


1.275 (11) 


-2.238(63) 


7 


1.1583(63) 


-2.263(43) 


1.161(17) 


-2.25 (12) 


1.159 (10) 


-2.259(70) 


8 


1.0620(58) 


-2.282(47) 


1.064(15) 


-2.27 (12) 


1.0626(95) 


-2.278(76) 


9 


0.9809(53) 


-2.299(50) 


0.982(13) 


-2.29 (12) 


0.9815(87) 


-2.294(83) 


10 


0.9117(49) 


-2.315(54) 


0.913(11) 


-2.30 (12) 


0.9122(81) 


-2.309(89) 



Table 3: Results for In(ALmax) with different fit strategies. 



of asF ~ 0.28. For A'^f = 2 a similar effect was visible only for larger coupling [11]^. 
These findings underline the necessity of going to weak coupling before applying 
perturbation theory. 

Clearly the present work has brought us a good step closer to the computation 
of the A-parameter in 4-flavor QCD, which may then be perturbatively connected 
to e.g. the 5-flavor MS coupling at the Z-pole. However, the technically introduced 
scale Ljnax remains to be expressed in physical units through large volume 4-flavor 
simulations. Apart from the challenge of tuning more parameters, one needs 
to treat a massive charm quark at small enough lattice spacing. Presently this 
appears to be a considerable challenge due to a severe slowing down of lattice 
simulation algorithms at small lattice spacings [35, 36]. 
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Tor asF ~ 0.45 a similar deviation is visible but there are no non-perturbative data points 
in between agF ~ 0.28 and asF ~ 0.45 to see better where this sets in. 
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Figure 4: The running coupling in the Schrodinger functional scheme. The gray 
vertical dashed line is only for the guidance of the eyes to show what the pertur- 
bation theory predicts for asF at our lowest energy. 

Appendix: The MC data 

We here list the results of our MC simulations. Each row refers to a separate 
simulation of about 50000 trajectories unless otherwise noted. Some of these 
simulations consist of independent replica (between 1 and 16). Measurements 
were taken after every trajectory for which we chose trajectory length r = 1. 
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1.3844(42) 
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1.1119 (10) 


0.00070(15) 


1.2074(32) 


0.012435(40) 






8.6 


0.131860 


0.99575(77) 


-0.00021(14) 


1.0678(25) 


0.010382(35) 






9.2 


0.131309 


0.90315(54) 


0.00094(11) 


0.9662(20) 


0.010176(33) 








Table 4: The raw 


data for L/a = 


4 and L / a = 


8. 
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2.3507(92) 
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2 904 r29'l 
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0.136443 
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0.00053 (11) 


2 529 f23) 


-0.000031(37) 


6.45 
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1.6948(46) 




-0.000294(94) 


1.953 (14) 


-0.000922(31) 
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0.133261 


1.2500(24) 
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0.132538 
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0.000435(62) 


1.2051(50) 


-0.000347(20) 
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0.131935 
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0.000154(57) 








9.45 


0.131411 


0.8975(12) 


0.000237(51) 


0.9628(31) 


-0.000547(16) 



Table 5: The raw data for L/a = 6 and L/a = 12. 
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0.000233(40) 
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-0.000666(15) 
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0.132736 
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0.000069(38) 
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0.000328(36) 
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-0.000552(12) 
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0.131442 
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0.000284(32) 


0.9961(48) 


-0.000504(11) 



Table 6: The raw data for L/a = 8 and L/a = 16. The run L/a = 16, 13 = 5.44 
has 98000 trajectories. 



19 



References 



[1] F. Butler, H. Chen, J. Sexton, A. Vaccarino, and D. Weingarten. Hadron 
masses from the valence approximation to lattice QCD. Nucl. Phys., 
B430:179, 1994, hep-lat/9405003. 

[2] S. Aoki et al. Quenched light hadron spectrum. Phys. Rev. Lett., 84:238, 2000, 
hep-lat/9904012. 

[3] A. Ali Khan et al. Light hadron spectroscopy with two flavors of dynamical 
quarks on the lattice. Phys. Rev., D65:054505, 2002, hep-lat/0105015. 

[4] S. Aoki et al. 2+1 Flavor Lattice QCD toward the Physical Point. Phys. 
Rev., D79:034503, 2009, arXiv:0807.1661. 

[5] C. Alexandrou et al. Light baryon masses with dynamical twisted mass 
fermions. Phys. Rev., D78:014509, 2008, arXiv:0803.3190. 

[6] S. Diirr et al. Ab-initio Determination of Light Hadron Masses. Science, 
322:1224, 2008, arXiv:0906.3599. 

[7] Martin Liischer. Advanced lattice QCD. Les Houches 1997, Probing the 
standard model of particle interactions, Pt. 2, 1998, hep-lat/9802029. 

[8] Martin Liischer, Rainer Sommer, Peter Weisz, and Ulli Wolff. A Precise 
determination of the running coupling in the SU(3) Yang-Mills theory. Nucl. 
Phys., B413:481, 1994, hep-lat/9309005. 

[9] Stefano Capitani, Martin Liischer, Rainer Sommer, and Hartmut Wittig. 
Non-perturbative quark mass renormalization in quenched lattice QCD. Nucl. 
Phys., B544:669, 1999, hep-lat/9810063. 

[10] S. Takeda et al. A scaling study of the step scaling function in SU(3) gauge 
theory with improved gauge actions. Phys. Rev., D70:074510, 2004, hep- 
lat/0408010. 

[11] Michele Delia Morte et al. Computation of the strong coupling in QCD with 
two dynamical flavours. Nucl. Phys., B713:378, 2005, hep-lat/0411025. 

[12] S. Aoki et al. Precise determination of the strong coupling constant in N(f) = 
2+1 lattice QCD with the Schrodinger functional scheme. JHEP, 0910:053, 
2009, arXiv:0906.3906. 



20 



[13] Thomas Appelquist, George T. Fleming, and Ethan T. NeiL Lattice Study of 
Conformal Behavior in SU(3) Yang-Mills Theories. Phys. Rev., D79:076010, 
2009, arXiv:0901.3766. 

[14] Q. Mason et al. Accurate determinations of from realistic lattice QCD. 
Phys. Rev. Lett, 95:052002, 2005, hep-lat/0503005. 

[15] M. Gockeler et al. A determination of the Lambda parameter from full lattice 
QCD. Phys. Rev., D73:014513, 2006, hep-ph/0502212. 

[16] E. Shintani et al. Lattice study of the vacuum polarization function and 
determination of the strong coupling constant. Phys. Rev., D79:074510, 2009, 
arXiv:0807.0556. 

[17] Fatih Tekin, Rainer Sommer, and UUi Wolff. Symanzik improvement of lat- 
tice QCD with four flavors of Wilson quarks. Phys. Lett., B683:75, 2010, 
arXiv:0911.4043. 

[18] Martin Liischer, Peter Weisz, and UUi Wolff. A Numerical method to compute 
the running coupling in asymptotically free theories. Nucl. Phys., B359:221, 
1991. 

[19] Particle Data Group. Particle Physics Booklet, http://pdg.lbl.gov, 2008. 

[20] Martin Liischer, Rajamani Narayanan, Peter Weisz, and UUi Wolff. The 
Schrodinger functional: A Renormalizable probe for non-Abelian gauge the- 
ories. Nucl. Phys., B384:168, 1992, hep-lat/9207009. 

[21] Stefan Sint. On the Schrodinger functional in QCD. Nucl. Phys., B421:135, 
1994, hep-lat/9312079. 

[22] Stefan Sint. One loop renormalization of the QCD Schrodinger functional. 
Nucl. Phys., B451:416, 1995, hep-lat/9504005. 

[23] Stefan Sint and Rainer Sommer. The Running coupling from the QCD 
Schrodinger functional: A One loop analysis. Nucl. Phys., B465:71, 1996, 
hep-lat/9508012. 

[24] Achim Bode, Peter Weisz, and Ulli Wolff. Two loop computation of the 
Schrodinger functional in lattice QCD. Nucl. Phys., B576:517, 2000, hep- 
lat/9911018. Erratum-ibid.B600:453,2001, Erratum-ibid.B608:481,2001. 



21 



[25] Rainer Sommer. Non-perturbative QCD: Renormalization, 0{a)- 
improvement and matching to heavy quark effective theory, hep-lat/0611020, 
2006, hep-lat/0611020. 

[26] Bernd Gehrmann, Juri Rolf, Stefan Kurth, and Ulli Wolff. Schrodinger 
functional at negative flavour number. Nucl. Phys., B612:3, 2001, hep- 
lat/0106025. 

[27] Achim Bode et al. First results on the running coupling in QCD with two 
massless flavors. Phys. Lett, B515:49, 2001, hep-lat/0105003. 

[28] M. Liischer and P. Weisz. 0{a) improvement of the axial current in lattice 
QCD to one-loop order of perturbation theory. Nucl. Phys., B479:429, 1996, 
hep-lat/9606016. 

[29] Giulia de Divitiis et al. Universality and the approach to the continuum limit 
in lattice gauge theory. Nucl. Phys., B437:447, 1995, hep-lat/9411017. 

[30] Stefan Sint. The Schrodinger functional with chirally rotated boundary con- 
ditions. PoS, LAT2005:235, 2006, hep-lat/0511034. 

[31] Stefan Sint. Lattice QCD with a chiral twist, hep-lat/0702008, 2007, hep- 
lat/0702008. 

[32] Jenifer Gonzalez Lopez, Karl Jansen, Dru B. Renner, and Andrea Shindler. 
Chirally rotated Schrodinger functional: non-perturbative tuning in the 
quenched approximation. PoS, LAT2009:199, 2009, arXiv:0910.3760. 

[33] Urs M. Heller. The Schrodinger functional running coupling with staggered 
fermions. Nucl.Phys., B504:435, 1997, hep-lat/9705012. 

[34] Paula Perez- Rubio and Stefan Sint. The SF running coupling with four 
flavours of staggered quarks. PoS, LAT2007:249, 2007, arXiv:0710.0583. 

[35] Luigi Del Debbio, Haralambos Panagopoulos, and Ettore Vicari. Theta de- 
pendence of SU(N) gauge theories. JHEP, 0208:044, 2002, hep-th/0204125. 

[36] Stefan Schaefer, Rainer Sommer, and Francesco Virotta. Investigating 
the critical slowing down of QCD simulations. PoS, LAT2009, 2009, 
arXiv:0910.1465. 



22 



